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Abstract 



We study the transient statistical properties of short and long Josephson junctions under the 
influence of thermal and correlated fluctuations. In particular, we investigate the lifetime of the 



■ superconductive metastable state flnding the presence of noise induced phenomena. For short 

O . Josephson junctions we investigate the lifetime as a function both of the frequency of the current 



driving signal and the noise intensity and we flnd how these noise-induced effects are modified by 



' the presence of a correlated noise source. For long Josephson junctions we integrate numerically the 

sine-Gordon equation calculating the lifetime as a function of the length of the junction both for 
' inhomogeneous and homogeneous bias current distributions. We obtain a nonmonotonic behavior 

of the lifetime as a function of the frequency of the current driving signal and the correlation time 
' of the noise. Moreover we find two maxima in the nonmonotonic behaviour of the mean escape 



time as a function of the correlated noise intensity. 



I. INTRODUCTION 



* e-mail: augello@gip . dft . imipa. it 

te-mail: valentid@gip.dft.umpa.it In recent years a great attention was paid 
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to the study of Josephson junctions (JJs) be- 



^Electronic address: |http: / /gip.dft.unipalt] cause of their use both as superconducting 

quantum bits , IJ] and nanoscale su- 

perconducting quantum interference devices 
for detecting weak magnetic flux change j^. 



This widely used device is very sensitive to 
magnetic flux change and it is composed by 
two coupling high-Tc JJs in a superconduct- 
ing ring. JJs are good candidates to real- 
ize superconducting quantum bits (qubits) 
for quantum information processing. In par- 
ticular, they were studied at very low tem- 
perature in devices makinguse of charge Q], 
flux [7| and phase qubits [8|]. In JJs, work- 
ing both at high and low temperatures, the 
environment affects strongly the behavior of 
the system. In high temperature supercon- 
ductors (HTSs) the presence of low-frequency 
noise, whose intensity is related to the fluc- 
tuations in the bias current, temperature and 

CI 

magnetic field, was experimentally found [9|. 
Also in the low temperature superconductive 
devices it is very difficult to avoid the influ- 
ence of environment that constitutes mainly 
a decoherence source for the system. In par- 
ticular, the effects on the coherence time 
of weak current noise in Josephson vortex 
qubits (JVQs), composed by an uniform long 
JJ, were investigated 0, 10]. 

In this framework the study of transient 
dynamics of JJs in the presence of noise 
sources is very interesting for the understand- 
ing of the interaction between these systems 
and environment. In particular, the effects 
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of noise strongly influence the current- voltage 
characteristic of JJs [ll,^, The dynam- 
ics of a JJ is studied considering a fictitious noise. For given values of frequency of the 



Brownian particle moving in a washboard po- 
tential IJ] and the behavior of the current- 
voltage characteristic of a JJ is strictly re- 
lated to the lifetime of the superconductive 
metastable state of the particle. The decay 
of the particle from the metastable state, in 
fact, depends on the fluctuations of the volt- 
age across the junction. Recently noise in- 
duced effects were experimentally observed 
in underdamped Josephson junctions 
and the switching to resistive state of an an- 
nular Josephson junction due to thermal ac- 
tivation was analyzed |l7|. 

In the present work, we study the tran- 
sient dynamics of short overdamped and long 
JJs under the influence of fluctuating bias 
current and oscillating potential. We nu- 
merically calculate the lifetime of the super- 
conductive metastable state also called the 
mean switching time (MST) to the resis- 
tive state for short JJ (SJJ) and long JJ 
(LJJ). We demonstrate the presence of noise 
induced effects such as resonant activation 
20| a nd noise enhanced stability 
122]. We analyze both the ef- 
fects of thermal and correlated noise sources. 
In SJJ and LJJ we consider white noise, ac- 
counting for the thermal fluctuations, and 
correlated (colored) noise separately. More- 
over, in LJJ we present an analysis consider- 
ing together the effects of white and colored 



(RA) yp 

(NES) y 
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driving signal and suitable noise intensity, we 
find maxima of the lifetime of the supercon- 
ductive state. This is an interesting feature 
for the study of the coherence time of these 
devices. Our results hold for low tempera- 
ture superconducting devices when we con- 
sider only the effects of the colored noise, and 
they can be extended to high temperatures, 
when both colored and white noise come into 
play. 

II. SHORT OVERDAMPED JUNC- 
TIONS 

A. Model 

The study of SJJs is performed in the 
framework of the resistively shunted junction 
(RSJ) model formalism [23]. To take into ac- 
count the transient dynamics of the system, 
a fluctuating current term in the RSJ model 
equation is considered. We obtain the follow- 
ing Langevin equation 24] 



with h the Plank constant. In Eq. ([T]) the 
time is normalized to the inverse of the char- 
acteristic frequency of the junction Uc- In our 
analysis ((t) is a colored noise source, gener- 
ated using an Ornstein-Uhlenbeck (OU) pro- 



cess 



25|, characterized by a correlation time 



Tc. The potential profile U((f)) of Eq. ([T]) is 
given by 



= 1 - COS0 - i{t)(j), 



(2) 



del) dU{(f)) 

-n = -^c—7- ^cCW, (1) 

at d(p 
where is the order parameter of the system, 
that is the phase difference of the wave func- 
tions in the ground state between left and 
right superconductive sides of the junction. 
The characteristic frequency of the Joseph- 
son junction is uJc=2eRNlc/f>', where e is the 
electron charge, R]^^ is the normal conduc- 
tivity, Jc is the critical current and h=h/2Ti 



where i{t) = io + f{t), iq = ib/Ic is 
the constant dimensionless bias current and 
f(t)=AsinLot is the driving current with di- 
mensionless amplitude A = is/Ic and fre- 
quency u {ih and ig represent the bias cur- 
rent and the driving current amplitude re- 
spectively). 

The dynamics of a JJ described by Eq. ([1]) 
is equivalent to the dynamics of a particle of 
coordinate moving in the washboard poten- 
tial U((f)). The motion is that of a Brownian 
particle because of the presence of the noise 
term [2^ . The OU process of Eq. ([1]) is rep- 
resented by the stochastic differential equa- 
tion l25i 



dm 



■-C{t)dt+ ^dW{t) 

Tr Tr 



(3) 



where 7 is the noise intensity and W(t) 
is the Wiener process with the usual sta- 
tistical properties: {dW{t)) = 0, and 
{dW{t)dW{t')) =6{t-t')dt. The correlation 



function of the OU process is ence of white noise (see also Ref. 24|). 



, . I In Fig. 1 we report the behavior of MST 
{C{t)C{t )) = T^e . (4) 

''c vs Lu, for two values of the noise intensity, 

We investigate the dynamics of the Brownian namely 7 = 0.02 and 7 = 0.5, and differ- 

particle in the presence of a time- dependent ent values of Tc. The non-monotonic behav- 

nonlinear periodic potential, by solving nu- ior of the curves shows that the RA phe- 

merically Eq. ([T]), for ujc=^- nomenon, already found in the presence of 



To study the lifetime of the supercon- white noise |20|, l2J| , appears also with colored 

ductive state we take, as initial condition noise. The values of MST around the mini- 

for the particle, the minimum of the po- mum are influenced by the variation of Tc, 

tential profile, corresponding to the condi- more strongly for higher values of the noise 

tion (f)o=arcsin(io) and we calculate the time intensity. Moreover, for higher intensity val- 

spent by the particle to reach the next max- ues (see Fig. [1]), we find, in a wide range of 

imum. We perform a number of simulations frequency (0.3 < u < 0.8), a non-monotonic 

ranging from = 5000 to = 10000 to esti- behavior of MST as a function of the cor- 

mate MST with a good approximation. The relation time. This is shown more clearly 

procedure is repeated for different values of in Fig. El where we report MST vs Tc for 

the system parameters such as the amplitude 00 = 0.7. In Fig. [3] we report the curves of 

A and frequency u of the driving current sig- MST and SD vs cu. We note a range of fre- 

nal, the correlation time Tc and the intensity quency (0.2 <uj < 0.8) in which MST and SD 

7 of the colored noise. present a minimum, which indicates the sup- 
pression of timing error effects, already found 

-r, T -n n in the presence of white noise l20i. Moreover 

B. Liietime or metastable state ' ' 

when we consider the RA phenomenon in 
Our analysis of the lifetime of the super- the presence of colored noise, we find a seal- 
conducting state concerns the behavior of the ing effect depending on the values of the 
curves representing MST and the correspond- correlation time. In Eq. (jl]) if the correla- 
ing standard deviation (SD) vs u, and MST tion time is greater than the characteristic 
vs 7, for different values of Tc. In the fol- time scale of the system, |t-t'|, we obtain 
lowing figures we present the curves of MST (C(t)C(t')) ^1 color ed/^Tc. By comparison with 
calculated in the presence of colored noise, in- the correlation function of the white noise, 
eluding also the curves obtained in the pres- {^{t)C,{t'))='ji^h.ite^{t — t'), we can define the 
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FIG. 1: (a): MST vs u for white noise and dif- 
ferent r^, 7=0.02. (b): MSTvs a; for white noise yIG. 2: (a): MST vs for 7=0.02. (b): MST 
and different r^, 7=0.5. In both panels io=0.8 ^^q.S. In both panels w=0.7, io=0.8 

and A=0.7. 



and A=0.7. 



effective intensity of the colored noise scaled 
by a factor 1/2tc and equivalent, in this ap- 
proximation, to the intensity of the white 
noise 

'ycolored 



(5) 



We check this equivalence by considering 
the effect on the system of a white noise with 
intensity ^white=^-^^'^ and a colored noise 
with intensity 7co«ored=0.5, setting rc=50. In 
Fig. H^, we report the curves of MST vs 



ijj calculated for these two values of noise 
intensity. We see that for MST less than 
50 there is a good agreement between the 
curves. In these conditions the evolution oc- 
curs with time scale smaller than the corre- 
lation time and the system is not affected by 
the noise memory. In Fig. Hb, we show the 
behavior of MST function of the noise 
intensity for different values of the correla- 
tion time. We note the appearance of noise 
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FIG. 3: (a): MST vs uj for different r^, 7=0.02, 
io=0.8 and A=0.7. (b): Standard deviation 
(SD) vs UJ. 
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FIG. 4: (a): MST vs w with io=0.8, A=0.7: 
white noise (hne) and colored noise with rc=50 
(dots), (b): MST vs 7 for white noise and col- 
ored noise with different values of Tc, w=0.5, 
A=l, io=0.5. 



III. LONG JUNCTIONS 



A. Model 

enhanced stability (NES), a phenomenon al- 
ready found in short JJs in the presence of The investigation of LJJs is developed in 
white noise 20|, l2J]. Here we observe a the framework of the sine-Gordon model 23l . 
range of values of correlation time, namely 0]. The dynamics of a LJJ is described by 
0.01<rc<l, in which the curves present a the motion of a phase string in a one di- 
non-monotonic behavior. For rc=5,10, the mensional washboard potential. The tran- 
non-monotonic behavior disappears. sient dynamics is represented by a nonlinear 



partial differential equation with a stochastic a fluctuating term given by the sum of a cor- 



term 
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28| related noisy current signal, icn{x,t), and a 



(9^0 d(f) d'^cj) thermal current signal, iwn{x,t), namely 

/5^+^+^ = ^(a^)-sin0+i/(a:,t), (6) 

if{x,t) =icnix,t) + i^ri{x,t). (8) 

where f3 = UcRC, with R and C the LJJ 

equivalent resistance and capacitance, re- The integration of Eq. ([6]) is performed us 
spectively. Here i{x) and if{x,t) are the m an imphcit finite-difference method 



bias current density and the fluctuating cur- 30], that allows to obtain a system of equa- 

rent density normalized to the critical cur- tions, whose solution is calculated by a tridi- 

rent density h of the junction, respectively, agonal algorithm. The integration step both 

We treat both the homogeneous bias cur- for the time and the space is 0.05. The num- 

rent density case, namely i{x) = h, and ber of simulation ranges from N = 100 to 

the inhomogeneous one, namely i{x) = ^ = 5000. 



{ibL)/{7ryJx{L — x)), where L = l/Xj with / 

the length of the junction and A j the Joseph- B. Effects of correlated noise 



son penetration depth [27|, 128|. In Eq. 
time is normalized to the inverse of the char- 
acteristic frequency, Uc = 2eRIc/h, of the 
junction. Analogously length is normalized 
to the Josephson penetration depth. The 
boundary conditions for Eq. are: 



Here we study the effects of a correlated 
noise on the transient dynamics of a LJJ, 
neglecting, in Eq. IQ, the term t) 
and considering a fluctuating current signal 
icn{x,t) with correlation function 



d<f){0,t) _ d<P{L,t) _^ {tUx,t)tUx',t')) = ^^^6{x-x')e-^, 
dx dx 

As initial position of the phase string in the where 2'ycn is the intensity of the correlated 

washboard potential, we consider the mini- noise. Therefore, the current density icn (2;, ^) 

mum, represented by the condition 0) = is subjected to fluctuations that are corre- 

arcsin{iij). lated in time (colored noise), showing a 5- 

First we present our analysis concerning correlated behavior in space (white noise) . In 

the effects of a fluctuating current signal de- Fig. we report the curves of MST vs the 

scribed by a correlated noise. Then we con- length of the junction, in the homogeneous 



sider the effects of both colored and white bias current case, for different va 
noise on the system, by inserting in Eq. ([6]) including white noise (tc = 0) 



23, 



ues of Tc, 



28|. We 



find a decrease of MST as tlie intensity of the 
colored noise •jcn increases (see Fig. ^ (7c„ = 
0.3) and Fig. (7cn = 0.7)). Moreover, we 
observe that, for fixed values of the noise in- 
tensity, MST increases as the correlation time 
becomes larger. We also present the numeri- 



in correspondence of L = 5. In order to in- 
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FIG. 5: (a): MST vs L with if,=0.7, 7cn=0.3. 
(b): MST vs L with ib=0.7, 7cn=0.7. 

cal results for the case of inhomogeneous bias 
current. In Fig. we show the behavior 
of MST function of the dimensionless 
length L. We observe that the curves present 
a non-monotonic behavior with a maximum 
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FIG. 6: (a): MST vs L with ife=0.7, 7c„=0.4. 
(b): MST vs oj with 7cn=0.4, A=0.7 and ^^=0.9, 

vestigate the presence of RA phenomenon, in 
Eq. we replace the term i{x) with a ho- 
mogeneous oscillating driving current signal 
given by i{t) = % + Asin{ujt), with A and 
LU amplitude and angular frequency, respec- 
tively, of the periodical driving signal. The 
curves obtained by numerical simulation are 
presented in Fig. and Fig. [7^. In both fig- 
ures we note the presence of a minimum in 
the curves of MST vs u. This is the signa- 
ture of the resonant activation phenomenon. 
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Moreover in Fig. [7^, in the frequency range and in the superconductive state in the pres- 
ence of a driving current signal (Fig.[8)D). The 
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FIG. 7: (a): MST vs u; with 7cn=0.9, A=0.7 and 
ife=0.9. (b): MST vs Tc with 7c„=0.9, A=0.7 and 
^6=0.9. 

0.3 < uj < 1.0, we find that curves with dif- 
ferent Tc are overlapping. In order to investi- 
gate in detail the behavior in this frequency 
range we report, for the same parameter val- 



curves of Fig. [HK present a maximum for dif- 
ferent values of Tc. In Fig. [8b we find, in all 
curves, a double peak. The comparison of the 
two panels in Fig. [8] shows that the driving 
signal causes the increase of the two maxima 
(see Fig. [8)d) already present in the absence 
of periodical term (see Fig. [Hti). The non- 
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ues of Fig. the MST as a function of r, ^^y_ ^^^.^^.^^ ^^^^^ 

for u; = 0.9, finding a nonmonotonic behav- 



ior. We also investigate the presence of NES 
for the homogeneous current case in the re- 
sistive state (Fig. [8^) without driving signal. 



{ib = 1.2) with L=0.5 (b): MST vs 7c„ in the 
superconductive state = 0.9) with L=0.5, 
UJ = 0.9 and A=0.7. 



monotonic behavior of MST function of 
the noise intensity 7cn indicates the presence 
of a NES effect. The curves of Fig. [S] have 
been calculated for L = 0.5. In Fig. [Hk, we 
report the curves of MST vs 7cn for L = 5, 
obtained in the presence of periodical driv- 
ing signal for the superconductive state. We 
note that, using the same noise intensities, 
for longer junction the decay time is shorter. 
This can be explained recalling that the ef- 
fect of the spatial correlation disappears on 
a longer distance. This corresponds to the 
presence of white noise, which is responsible 
for faster dynamics and, then, shorter escape 
time. 

C. Effects both of correlated and 
thermal noise 

In this section we investigate the effects 
both of thermal and correlated noise, con- 
sidering in Eq. the fluctuating current 
term if{x,t) given by Eq. (jS]). In Fig. [Hb, we 
present the curves of MST vs 7c„ for different 
white noise intensities 7^„, with L = 5. 

We note that when is greater than 
7cn the effects of colored noise disappear. If 
we consider 7^^ suitably lower than 7c„, the 
effects of colored noise become evident. In 
fact, when 7c„ is at least five order of mag- 
nitude higher than 7^„, the curve of MST 
vs 7cn (empty triangles) matches with good 
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FIG. 9: (a): MST vs 7cn in the superconduc- 
tive state {ih = 0.9) with L=5, for different val- 
ues of correlation time, (b): MST vs 7cn in the 
superconductive state {ib = 0.9) with L=5, for 
different values of white noise intensity. 

agreement that (black circles) obtained in 
the absence of white noise (see Fig. [9]d for 

7^n = 10-«). 

IV. CONCLUSIONS 

We studied the effects of white and colored 
noise on the transient dynamics of short and 
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long Josephson junctions. We investigated dom fluctuations influence the mean switch- 

the lifetime of the metastable state finding ing time of short and long Josephson junc- 

noise induced effects, namely resonant activa- tions. We found that these times are strictly 

tion and noise enhanced stability. The results affected by the characteristic parameters of 

obtained throw light upon the role played by the system such as intensity of noise, fre- 

different noise sources in the dynamics of su- quency of the driving signal and correlation 

perconductive devices, explaining how ran- time of the colored noise. 
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